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Abstract Based on the method which is given in Ref. [Sun et.al. arXiv:0904.0092vl], we 
present another 9x9 unitary ^—matrix, solution of the Yang-Baxter Equation, is obtained in 
this paper. The entanglement properties of matrix is investigated, and the arbitrary degree 
of entanglement for two-qutrit entangled states can be generated via ^-matrix acting on 
the standard basis. A Yang-Baxter Hamiltonian can be constructed from unitary ^—matrix. 
Then the geometric properties of this system is studied. The results showed that the Berry 
phase of this system can be represented under the framework of SU(2) algebra. 
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1 introduction 



Quantum entanglement(QE), the most surprising nonclassical property of quantum systems, 
plays a key role in quantum information and quantum computation processing! 2 3 , 4 5 ]. Be- 
cause of these applications, QE has become one of the most fascinating topics in quantum 
information and quantum computation. On the other hand, the geometrical phase |'6], such 
as Berry phase(BP), is another important concept in quantum mechanics 17,8,9 JTOlfTTl . In 
recent years, a lot of works have been attributed to BP[ 12 1, because of its possible appli- 
cations to quantum computation(the so-called geometric quantum computation) 1 1 3 II flll 151 . 
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Such concern is motivated by the belief that geometric quantum gates should exhibit an in- 
trinsic fault tolerance in the presence of some kind of external noise due to the geometric 
nature of the BP. 

Yang-Baxter Equation(YBE) II 1 6 II 1 7lfT8l was originated in solving quantum integrable 
models, but recently has been shown to have a deep connection with topological quantum 
computation and entanglement swapping fT9ll20ll21ll22ll23ll24ll25ll26l . In Ref.|27) , the au- 
thors point out YBE can be tested in terms of quantum optics. In a very recently work[26|, 
it is found that any pure two-qudit entangled state can be achieved by a universal Yang- 
Baxter Matrix assisted by local unitary transformations. However, the solution R(x)— matrix 
in Ref.| 26 ] only dependent on one parameter. So we can't construct a Yang-Baxter Hamilto- 
nian as in Ref.[24 28 1. In this paper, we obtain a time-dependent solution of YBE, R(x, <pi , 92) 
<P\ and 92 are time-dependent, so we can construct Yang-Baxter Hamiltonian. Consequently, 
we can study entanglement properties and Berry phase for this system. 

This paper is organized as follows: In Sec[2] we present a 9 x 9 Yang-Baxter matrix. 
By means of negativity, we investigated the entanglement properties of R(x, <p\ , 92)-matrix. 
We show that the arbitrary degree of entanglement for two-qutrit entangled states can be 
generated via the unitary matrix R(9, q>\, 92) -matrix acting on the standard basis. In Secj3] 
we construct a Hamiltonian from the unitary ^(0, <pi , <p2)-matrix. The Berry phase of the 
system is investigated, and the results showed that the Berry phase of this system can be 
represented under the framework of SU(2) algebra. The summary is made in the last section. 

2 Unitary solution of Yang-Baxter Equation and its entanglement properties 

The usual YBE takes the form, 



The spectral parameters x and y which are related with the one-dimensional momentum 
play an important role in some typical models| 16 1. The asymptotic behavior of R(x, <p\ , 92) 
is x-independent, i.e. limRij + \ (x, 91, 92) <* where are braiding operators, which satisfy 
the braiding relations, 



where the notation b t = is used, &;,,•+] represents lj ® I2 ® I3 • • • ® Si.i+i ® * ' * ® In > 
and \j is the unit matrix of the j-th particle. 

As is known, Hecke algebras are intimately connected with braiding groups. In fact, 
braid algebra is subalgebra of Hecke algebra. And we can construct a representation of 
braid algebra from Hecke algebra. A unitary solution of YBE can also be constructed 
from a representation of Hecke algebra. Let us review Yang-Baxterizationl 29 . 30 11281 of 
Hecke algebra. M,, a Hermitian matrix(/.e.M/ = M,), satisfies the Hecke algebraic relations: 
MiMi+iMi + gMi = M i+ \MiM i+ i + gM i+ \ and Mf = aM t + jS/;. For convenience, we set 
a = 1 and /3 = g = 2. let the unitary Yang-Baxter matrix take the form, 



Substituting Eq® into EqQ, one has F(x) + F(y) + F(x)F(y) = [1 + 2F(x)F(y)]F(xy). 
The unitary condition (i. e. , Rj (x) = RJ~ 1 (x) = Ri (x~ 1 ) ) can be tenable only on condition that 



Ri(x)R i+ i {xy)Ri{y) = R w {y)Ri{xy)R i+1 (x) 



(1) 




(2) 



R i (x)=p(x)[li + F(x)M i \ 



(3) 



3 



F(x) +F(x- V ) +F(x)F(x- 1 ) = and p(x)p( X - l )[l + 2F(x)F( X - 1 )} = 0. In addition, the 
initial condition Ri(x = 1) = 7, yields F(x=l)=0 and p(x = 1) = 1. Taking account into these 
conditions, we obtain a set solutions of F(x) and p(x), 

In this paper, we choose basis {| 11), 1 10), |01), 1 1 - 1),|00),| - 11), |0 — 1),| - 10), | - 1 - 
1)} as the standard basis. Based on calculation, a 9 x 9 matrix M which satisfies the Hecke 
algebraic relations is realized as, 

M cd = \a^d + q2Sabo\a^d 

+ Q~ ' S a b- 1 ajicytd + 9 1 1 S cc ll c^a+b 

(4) 

+ ?2 1 ^cdO I c^a^b + QScd-l I cjtayib 
+ SadSbc \a^b 

Where q\ = e"?' , q 2 = e"? 2 and Q = qiq 2 , with the parameters (p\ and <p 2 both are real. The 
denotes M"^ = M a /, cl j are used. The denote 8 a bc = 1, if and only if a = b = c; otherwise, 
the denote 8 a b c = 0. This solution is not equivalent to the solution in Ref. |28l . Substituting 
Eq([4ji into EqlfJ), the unitary solution of YBE can be obtained as following, 

R(x, <Pu92)td = P (x) [Sabcd +F(x)M%] (5) 
The matrix form of Ri(x, <Pi , 92) can be recast as, 
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where &=x 1 — x, b — 2x+x 1 . The Gell-Mann matrices, a basis for the Lie algebra SU(3)| 3 1 1, 
% u satisfy [7^,/^] = i'/^^ v 7 v (A,^l, V = 1, ■ ■ -,8), where 7^ = For the later convenience, 
we denote 7^ by, 7± = I\ ± /Y2, V± = V4 =f iV5,E/± = 76 ± 1/7, F = -^7g. In this work, we get 
rise to three sets of realization of SU (3) as: 

t/i'Wf^, v^W+t/f, 

jW= 1 (F 1+ y 2 )_ 2/3/3 _iy l72; 
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'i { l ] = u±u*, u { i ] =v±%, vi 2) =ifv^, 
4 2) = \[-Wl-ii)+k ( Y i - Y 2) +iiY 2 - Y X I\] , 

7 (2) = _ [ 1 (/ 3 + ,3) + 1 (7] + y 2 ) + 2 /3/ 3 + 1 y lYl] . 



4 3) = ll-Wl-fy - l(Yi-Y 2 ) +I\Y 2 - ¥,% 
{ 7(3) = !(/?+/!) - i(F! +F 2 ) - |/?/| - lF,F 2 . 

We denote /« = /<*> ± tff , = vf > T iV s ( Vj? = ^ ± ^ • 7 W = ^4*' (* = 1, 2, 3). 

These realizations satisfy the commutation relation [ffl ,ljP] = iSijfxuvIv (^)Mi v = ^ ' ' 
• , 8; /, y = 1 , 2, 3) . So the whole tensor space C 3 <8> C 3 is completely decomposed. In addition, 
each block of ^-matrix can be represented by fundamental representation of SU(3) algebra, 
i.e. C 3 ®C 3 = C 3 ®C 3 ®C 3 . 

For i-th and (i + l)-th lattices, ^-matrix can be expressed in terms of above operators, 

R(0, <p u (p 2 ) = ±a[I^ ] +/ m + Q{V W + U { + 1) ) 

+ G -i {c/ n) +v o) )+/ (2) +/ (2) 

+I ( 3 ) + p) +q2{v ^ +u (jy } 
+ 92 - 1 (y! 3) + f/| 3) )] + f(/®7). 

We can introduce a new variable with x=e' e , and 6 may be related with entanglement 
degree. When one acts R(9,<p\,(p2) on the separable state \mn) , he yields the follow- 
ing family of states \y) mn = £j^!/&w|m/t)(m,n=l,0,-l)- For example, if m=l and n=l, 
|^)n = i(/7|ll)+a ? - 1 |0-l)+a^ 1 |-10)).By means of negativity] 32 , 33 34], we study 
these entangled states. The negativity for two qutrits is given by, 



N(p)i 



llp^ll-l 



(7) 



where ||p rA || denotes the trace norm of p TA , and p Tk denotes the partial transpose of the 
bipartite state p. i.e., {p) l f \ = (P TA )i A h- ^ n ^ act ' ^(P) corresponds to the absolute value of 
the sum of negative eigenvalues of p Tk , and negativity vanishes for unentangled states 1331 . 
Then we can obtain the negativity of the state 1 1 as 



N(6) = - (s in 2 6 + | sin 6 \ \/ 1 + 8cos 2 6 ) . 



(8) 



When \a\ = \b\, namely x = e'?, the state \xj/)u becomes the maximally entangled state 



of two qutrits as \\jf)i\ = — lq x 1 1 — 1) — iq l 1 — 10)). In general, if one acts 
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the unitary Yang-Baxter matrix R(x) on the basis {|11), |10), |01), |1 - 1), |00), | - 11), |0- 
1), | — 10), | — 1 — 1)}, he will obtain the same negativity as Eq®. It is easy to check that 
the negativity ranges from to 1 when the parameter 9 runs from to n. But for 9 6 
[0, tc], the negativity is not a monotonic function of 9. And when x = e'7 , he will generate 
nine complete and orthogonal maximally entangled states for two qutrits. The QE doesn't 
dependent on the parameters (pi and q>2- So one can verify that parameter (pi and 92 ma y 
be absorbed into a local operation. Base on numerical calculation, the universality of YBE 
is proved by Chen et.al. in Ref . 1 26 1 . This unitary solution of YBE can generate entangled 
states, this solution may be a universal quantum gate. 



3 Yang-Baxter Hamiltonian and BP 

A Hamiltonian of the Yang-Baxter system can be constructed from the R(9, (p\ , 92)-matrix. 
As shown in Ref. [24], the Hamiltonian is obtained through the Schrodinger evolution of 
the entangled states. Let the parameters 9, be time-dependent as <p, = co,f. The Hamiltonian 
reads, 

at 

= 0# W , (9) 

k=\ 

where the superscript k denotes the k-Xh subsystem. The &-th subsystem's Hamiltonian ffw 
can be obtained as following, 

//(i) = c(l)[^™e(/| 1) + ^sin9Y^ 

(10) 

_ ^* e( yW + [/«) + ^ ibQ -i(yW +E /W)] 
= C (2)[-^s/«e(/l 2) +/i 2) ) - ^smBY®] 

(11) 

ff(3) = C(3)[-^*to0(4 3) +/i 3) ) " ^sin9Y^ 

(12) 

+ $iPq?<JJ® +yl 3 )) _ $i^(vf + tf< 3 >)] 

Where C(l) = - 4 ^'» e , C (2) = - 4 ^' "" e , C(3) = - 4 ^'^"" e and Q = ah + 0*. In 
terms of/?' (A = 1,2, • • ■ , &;k = 1,2,3), the Hamiltonian can be recast as following, 

g 

ffW=C(fc)£flf/f. (13) 
1=1 
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Compare Eq{R)]l, EqO, Eq{T2]i with Eq(l3l. one can obtain B^ as following 



= -^-sin9 cos (o(l)t + ^cos9smco(l)t 
= ^sin9smco(l)t + ^cos9 cos co(l)t 
= -^sind cos co{l)t + ^cosdsmco(l)t 
= ^sm0sinto(l)f + ^cos9 cos d)(l)t 



^sind ; B® = B® =0 



^ sinO cos co{i)t + ^cosO sin 0)(i)t 
^sin9 sin w(i)t- &cos9 cos 0>(z)f 



= ^s!n0costa(i)f+ ^cos9 sin o)(i')f 



■^j , in0sina)(i)f — ^-cosB cos co(i)t 



_ </2 



sinO 



Where i=2,3- The denotes 0)(1) = 12, a)(2) = Oi and 0)(3) = 0)2 are used. The Hamilto- 
nian for the k-th subsystem, ff(B(f)(*))W, depends on the parameters M (A = 1,2, • •• , 8;), 
which are the components of a vector B^. And B^ are a set of time-varying parameters 
controlling the k-th subsystem's Hamiltonian. After time f", Hamiltonian returns to its 
original form, i.e. #(B(0)) w = H(B(T^)) . According to this, one can easily verify peri- 
ods of the subsystems are rW = 2n/n, r< 2 ' = 27r/a)i and = 2n/o>2. The eigenstates 
of the first subsystem are found to be, 

\E^) =iv| 1) [/| 1) (| 10) + |01>) +e- in '\ - 1 - 1)] 

I^H^Hioj + loi)) 

=jvW(/W(|io> + 101)) + e - ;i2 '| - 1 - 1)). 

with the corresponding eigenvalues = sin 9, Eq ^ = and E_ = — sin 0. 

Where 



a4 1} = JMl^sine and 



^n$6y^ • ^ or me secon d and the third subsystems, the 
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eigenstates are found to be, 

1> + |-10»] 
1} + |-10»]. 

l)+l-ll))] 
l)+l-ll»]. 

with the corresponding eigenvalues E+ = ^^hco(i) sin0, E^ =0 andZji'' = — ^^hco(i) sin0. 
Where = y / 3±2 ^'" e and f£> = ^'"yjgf 1 ■ According to the definition of the BP (6), 

when the parameter [ s slowly changed around a circuit on the sphere of direction, then 
at the end of circuit, the eigenstates |£*)(a = +,0, — )evolves adiabatically from to T^ k \ 
the BP accumulated by the states \E„) are, 

^=if\E^\l\E^)dt. (14) 

Substitute these eigenstates into Eqdl4b. one can obtain the BP for these eigenstates(all 
phases are defined modulo 2k throughout this paper), 

' y+ = -{\-^rsind)2n 

< i = (15) 

. f- = {\-^rsind)2n 

In fact, the BP of this system can be represented under the framework of SU(2) alge- 
bra. First, we can introduce three sets SU(2) realizations in terms of three new sets of 
operators(Eqi[7]i,Eql|7]l,Eqi[7]l), 

' s f = ^(y« + [ /f) 

< S W = -±(vi k] + U {k] ) (16) 
They satisfy the algebraic relations of SU (2) group: 

[S { i\s U] ] =25 U S^, [S®,S { i ) ]=±8 ij S { i\ (5^)2 = 0(1,7 = 1,2,3), withS^ = s[ k) ± iS { 2 k] (k = 

1,2,3). By the way, their second-order Casimir operators are = \ (S+ S_ ' -\-S_ S+ ) + 

{sf 1 ) 1 . One can verify that the eigenvalues of J?^ are 5(5 + 1) = f and 0(0+ 1) = which 
correspond to spin- 1/2 system and spin-0 system. 



£ j 2 > ) = N { + ] [/| 2) 1 1 ) + e- ia » ' ( 1 - 

4 2) > = 73(-|0-l> + |-10» 

E m } = N m [/p) 1 1 1) + e -;<»if (|o - 

£ j 3) ) = /v| 3 ' [/j 3 > 1 00) +e-''^'(|l- 

4 3) ) = 73(-|l-l) + |-H)) 
£i 3) ) = /vi 3) [/1 3) 1 00) + e- iah ->{\ 1 - 
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When one substitutes these realizations into EqdlQb-Eq d 12b . he can recast Hamiltonian 
of the subsystems in terms of SU(2), 

H® = C{2)[\{B®sf +sf S^+B^S®} (17) 
ff( 3 ) = C(3)[i(BL 3) S« +B®S®)+Bf >S< 3 >]. 

Where = (B^)* = -^feV^B^ = ^sin6>,B W = (fl^)* = ±ib* e - im ^> andB^ = 
— 2^psin0(i=l;2). a)(i) and C(3) are defined below Eq(|13l>. So we can say the whole sys- 
tem equivalent to three spin-i subsystems and three spin-0 subsystems. In fact, we can 
introduce a time-independent 9x9 orthogonal matrix O (see Appendix A). By means of 
O, the whole system's Hamiltonian H and Casimir operators jfw are transformed into 
block-diagonal matrices, i.e. H = OHO T and = O ^^0 T are block-diagonal ma- 

trices, where T denotes the transpose of matrix O. For the subsystem 1, from Eq(|23l> we 
can get its Hamiltonian = H\' (BHq K For Ha~\ the eigenvalue of Casimir operator 

2 

is 0, and the BP is 0. So we can say the subsystem Hamiltonian Hq is equivalent 
to a spin-0 subsystem. For H ] , one can introduce two transformations, cosa = ^lAsinQ 
and cosB = -™ecosnt+3c s8sinS2, with a £ ( arccos h£L a rccos- ^) and B e [0,24 a is 

\/9-8sin 2 e i 

time-independent, and B is time-dependent. By means of this transformation, the Hamilto- 
nian H\' can be recast as H^p = C{l)(sinacosBS\ + sinasinBSi + cosaSj,). We substitute 

2 2 

these transformations into Eq d 15b . the BP of the subsystem ' can be recast as, 

2 

Y P = ^n(l - cosa) = t^(C)/2 (18) 

where Q(C) = 2n(\ — cosa) is the familiar solid angle enclosed by the loop on the Bloch 
sphere, and the parameter a comes from 9 which comes from the Yang-Baxterization of the 
unitary braiding operator. So the BP dependent on spectral parameter. Under the new basis 
Eqd22l. the eigenstates \E±*) can be recast as following (we neglected the global phase 
factor), 

14')) = _ e -<75 sm ||i )+cos ||2), (19) 
ff,., ,R . a, 



\E W ) = cos ~ 1 1) + e'P sin " |2) , (20) 



where |1) = 0-W(| 10) + 1 01) ) and 2) = 0\ — 1 — 1). As is known to all, they are spin 



coherent states. By means of Eq(l24li. the states \E±) can be recast as following, 

\E { ?)=exp[Zsf-Z*S { -} ) ]\2), 

\E W ) =exp[0 ( l ) -CS { l ) ]\\), (21) 
Where t, = e~® a /2. BP for spin coherent states has been investigated in Ref[35 1. So we can 

2 



say the subsystem H\ is equivalent to a spin- 4 subsystem. By means of the same method 
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the Berry phases for subsystem 2 and 3 may be obtained, y± = T^(l — cosa) = ^Q. (C) /2 
and ^ = 0. The whole system is equivalent to three spin— 1 /2 subsystems and three spin- 
subsystems. This Yang-Baxter Hamiltonian system is equivalent to the Hamiltonian in 
Ref.il). 



4 Summary 

In this paper, we have presented a 9 x 9 M-matrix which satisfies the Hecke algebraic rela- 
tions and derived a unitary R(9, <pi, 92) -matrix via Yang-Baxterization of the M-matrix. In 
the following, we show that the arbitrary degree of entanglement for two-qutrit entangled 
states can be generated via the unitary R(8, q>\, 92) matrix acting on the standard basis. Then 
the evolution of the Yang-Baxter system is explored by constructing a Hamiltonian from the 
unitary R(d,(pi, 92)-matrix. In addition, the BP of the system is investigated. By means of 
decomposition of the tensor product, the Berry phase of the whole system is explained. The 
whole system is equivalent to three spin— | subsystems and three spin-0 subsystems. Berry 
phase of this system is represented under the framework of SU(2) algebra. 

Thermal entanglement in multi-body system is an interesting and nature type of QE, so 
it is a good challenge to study the thermal entanglement in multi-body system. Yang-Baxter 
Equation is an important tool in this domain. 
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A Block-Diagonalize H and ^ w 

The time-independent 9x9 orthogonal matrix O reads, 
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The orthogonal matrix O satisfies the relation OO t = O t O = Ig x g, where O t denotes the transpose of matrix 
O. 

The orthogonal matrix O transforms the standard basis {| 1 1), 1 10), |01), 1 1 - 1),|00),| - 1 1), |0 - 1), | - 
10) , I — 1 — 1 ) } into a new set of basis. The relations of new basis and old basis are, 

' |l) = O^(|10) + |01» 

< |2) = 0|-1-1) 

. |3) = O-±-(-|10) + |01)) 
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' |4) = 0|11) 

, |5> = O^(|0-l) + |-10» 
6> = 0^HO-1} + |-1O» 
|7) = 0|11) 
, |8) = 0^(|1-1> + |-11» 
|9)=0^(-|1-1) + |-11» 



(22) 



{| 1), |2), 1 3), |4), 1 5), |6), |7), 1 8), |9)} are a set of new basis. By means of this set basis, the Hamiltonian H can 
be recast as block-diagonally form, 



H = OHO' 



(23) 



k=l 

where S W = ffi , ' 's are 2 X 2 matrix, and ff^' are 1 x 1 matrix with ff ^ = (0) . Under the new 

basis, the (3 x 3)-dimension matrix is decomposed into six blocks. 
Three sets of SU(2) realizations )16t can be recast as, 



f (2i 



1){2|, S 0) 


= |2><1|, 


= j(|1><1|- 


-|2><2|); 


4) (5 1, SL 2) 


= |5){4|, 


= i(|4><4|- 


-|5><5|); 


7) (8 1, S* 2 ' 


= |8}<7|, 


= l(\W\ 


-|8)(8|). 



(24) 



(25) 



(26) 



The seconde-order Casimir operators are Jf^> = |(|1)(1| + |2)(2|), Jl 
!(|7>(7| + |8>{8|). 



(2) 



»(|4){4| + |5)(5|), / 



0) : 
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